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Motivated by recent angle-resolved photoemission spectroscopy (ARPES) measurement that su- 
perconducting YBa2Cu307_a (YBCO) exhibits a d x 2_ s 2+s-symmetry gap, we show possible quasi- 
one-dimensional modulations of local density of states in YBCO. These aniostropic gap and defect 
induced stripe structures are most conspicuous at higher biases and arise due to the nesting effect 
associated with a Fermi liquid. Observation of these spectra by scanning tunneling microscopy 
(STM) would unify the picture among STM, ARPES, and inelastic neutron scattering for YBCO. 
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PACS numbers: 74.70. Pq, 74.20. Rp, 74.25.Ld 

After more than one and half decades of the discovery 
of high-temperature superconductors, there lies still one 
central question: are these materials behave like a Fermi 
liquid? Recent scanning tunneling measurements (STM) 
however have provided a new platform in detecting this 
critical issue, to which modulations of local density of 
states (LDOS) are observed in superconducting state of 
Bi 2 Sr 2 CaCu 2 8 (BSCCO) [1-3]. More recently, experi- 
mentilists managed to obtain Fourier transform (FT) of 
these LDOS modulations, which then can be immediately 
used to test whether the system is Fermi-liquid like. 

Analogous to angular resolved photo-emission spec- 
troscopy (ARPES), FT-STM can study electronic prop- 
erties connected to entire Brillouin zone of the system. 
When a system behaves as a Fermi liquid, Fermi sur- 
face (FS) scattering should be clearly identified in nor- 
mal state. While in superconducting state, it should also 
give information on how the symmetry of Cooper pair or- 
der parameter develops on the FS. Assuming the effect of 
defect, experimentally observed STM LDOS modulation 
can be ascribed to quantum interference arising because 
quasiparticles are scattered by the defect locally. The oc- 
currence of peaks and their evolution as the bias change 
in the FT-STM spectra seem to be explained well by 
Fermi liquid models [4-6], whose FS topology is consis- 
tent with previous ARPES results. Moreover, the origin 
of peak moving in FT-STM is very similar to what of 
the incommensurate peaks observed in inelastic neutron 
scattering (INS) [7]. Apparently this scenario is in great 
contrast to the stripe picture [8] , interpreted in the STM 
data of Howald et al. [2] and the one-dimensional (ID) 
INS data of Mook et al. [9]. 

Perhaps of technical involvement, the STM data done 
for YBCO so far have been limited to the CuO chain lay- 
ers only [10]. It becomes naturally to ask what would be 
seen by STM on Cu0 2 plane layers of YBCO . As is well 
known, due to the coupling to the ID CuO chains, the 
Cu0 2 plane is orthorhombic rather than tetragonal. The 
perfect D^h symmetry is broken which accounts naturally 
for the additional symmetry mixed in the order param- 



eter of the Cu0 2 plane. In this regard, it is commonly 
accepted that there is a subdominant s-wave component 
in addition to the dominant (£,.2_„2-wave superconduct- 
ing gap for Cu0 2 planes. This view is supported by ear- 
lier tunneling data [11] and verified by recent data of 
ARPES [12]. The ID like incommensurate INS peaks [9] 
can also be interpreted by this picture associated with 
a Fermi liquid [13]. When fitted to a simple d x 2_ y 2+s 
model, ARPES data reveals that the s-wave component 
is about 20 % of d-wave component [12]. 

Here, based on a Fermi liquid model, we investigate 
the LDOS specific to Cu0 2 planes of superconducting 
YBCO. Similar to other theoretical works [4-6], we con- 
sider the effect of defect, while pay special attention to 
the gap anisotropy. Our goal is to bring out the most 
conspicuous features and inspire STM measurements to 
be performed. Successful observation of our predictions 
in either real-space or FT LDOS could unambiguously 
identify (1) the goodness of Fermi- liquid behavior and 
(2) the d x 2_ y 2 + s symmetry of gap for YBCO. Further- 
more this could lead to a unified picture among ARPES 
[12] , INS [9] , and STM on YBCO. 

To investigate the LDOS spectra, we start by consid- 
ering the following Hamiltonian 



H — Hqcs + Hi, 
where the usual BCS part is 



(1) 



(2) 



and the part associated with the presence of one single 
extended impurity located at site is 



Hi = Yl % c *W + [ 

+ V (4 T c 0T + 4| c oi) • 



H.c. 



(3) 



The parameter St corresponds to local deviation of hop- 
ping due to impurity , 5A corresponds to local devi- 



ation of pairing gap due to impurity, and Vo is spin- 
independent impurity potential. Any spin-dependent im- 
purity potential ignoring spin-flip will have no effect in 
the lowest-order approximation. For extensive but weak 
impurity scattering, we consider local deviations to those 
terms (St±, SAi) couple the impurity site and its nearest 
neighbors and terms (<5^2, 8A2) couple impurity's near- 
est neighbors and next nearest neighbors. Similar kinds 
of Hamiltonian have been successfully used by Tang and 
Flatte [4] to explain the resonant STM spectra for Ni im- 
purities in BSCCO and by Wang and Lee [5] and Zhang 
and Ting [6] to explain the energy-dependent modulation 
of FT-STM spectra on BSCCO. We shall comment how 
our results are modified and their experimental relevance 
when more than one defect are present later. 

The local density of states measured by STM is given 

by 



D(r,uj)= - — Im[Gn(r,r, w + i0 + ) 

-Gn(T,T,-U-iO+)], 



(4) 



where Gy is the element of the 2x2 equal- location single- 
particle Green's function matrix G in Nambu representa- 
tion. The full G can be calculated via the usual Gor'kov- 
Dyson equation 

G(r, r, lu) = G°(r, r, u ) + £ G°(r, r', cu)f(r', r",cu) 

r' ,r" 

xG°(r",r,c), (5) 

where G°(ri,Tj) = G°(ri — rj) correspond to the no- 
defect ones of G(i"i, r^) and the sums (and the T-matrix 
T(r', r", u>)) in Eq. (5) are over the sites associated with 
Hi in Eq. (3). For weak impurity scattering, we shall 
calculate G in Eq. (5) up to first order of dt's, 6 A's, and 
Vo, i.e., in the Born limit or first order T-matrix approx- 
imation [6] . In our calculation, we use a 800 x 800 square 
lattice sites and the impurity is located in the center. For 
simple but reasonable reason, we choose 2Sti — A8t2 — 
— 2(5Ai = — 4<5A2 = Vo and assume these scales are samll 
such that the first order T-matrix approximation is valid. 
It should be emphasized that the weak impurity scatter- 
ing limit is enough to bring out essential physics. Strong 
impurity is not expected to give qualitatively new feature 
away from its neighborhood of interest (see later). For £k, 
we use a tight-binding band £k = ti(cosk x + cosk y )/2 + 
t^coskxCosky + ts(cos2k x + cos2fc a )/2 + f4(cos2fc :c cosfcy + 
cosk x cos2k y ) /2+t5COs2k x cos2k y (lattice constant a = 1), 
with ti_ 5 = -0.60,0.16,-0.05,-0.11,0.05 eV respec- 
tively for YBCO. The fine structure of the band do not al- 
ter the main results reported in this paper. For supercon- 
ducting gap of YBCO, we assume Ak = Ad/2(cosk x — 
cosky) + A s with Ad = 36.7meV and A s /Ad = s = 0.2. 
This gives the maximum gaps at antinodes (tt, 0) and 
(0,tt), A^ ax = 44meV and A^ ax = 29.3meV. Thus 
A^ ax /A^ lax = 1.5 in accordance with ARPES measure- 



ment [12]. In our calculation, we also introduce a fi- 
nite lifetime broadening 7 = 2 meV to the quasiparticle 
Green's function to smooth our data points by replacing 
u> + i0 + with lu + i"f. 



(c) QJ = -15meV (d) QJ = ISmeV 





(e) &J = -23meV (f) qj = 23meV 




(g) QJ = -23meV (li) QJ = -23meV 




s = 0.1 



s = 0.05 



FIG. 1. Bias energy dependent FT-LDOS plotted for the 
whole first Brillouin zone for a d x 2_ y 2 + s-wave superconduc- 
tor. All frames are for anisotropy s = 0.2 except for (g) and 
(h) with s — 0.1 and 0.05, which are intended to be compared 
to (e). Wavevectors qi and in (e) are the nesting ones 
shown in Fig. 2. 

Theoretically D(q, to), the Fourier transform of D(r, to) 
is more readily obtained. In Fig. 1, we calculate the bias 
energy dependent D(q,u>) over the whole first Brillouin 
zone. Energies have been taken from — 23meV to 23meV 
— a scale less than the maximum gap at antinode (0, 7r) of 
^max = 29.3 meV. All frames are for anisotropy s = 0.2 



except (g) and (h) with s = 0.1 and 0.05, which are 
intended to be compared to (e) . The case of small s re- 
covers the results of Wang et al. [5] and Zhang et al. [6] 
for a pure d x 2_ y 2-wave superconductor. When energy is 
low [case (a) and (b)], only node-to-node scattering is al- 
lowed, to which the finite q interference is not prominent. 
The most important interference appears to those zero or 
small q scatterings, which nevertheless exhibit clear fea- 
ture of anisotropic gap. When energy is high [case (e) 
and (f)], in contrast, finite q interference peaks emerge 
in addition to the strong interference structure centered 
at q = 0. 
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FIG. 2. The Fermi surface (FS) used in our calculation 
(dash lines) and contours of excitation energy Ek = 23 meV 
(solid lines) are shown. Nesting wavevectors qi and q' x con- 
necting locally almost parallel contours are intimately related 
to interference peaks. While q 2 is shown as an example of no 
assocaited interference peak. For better visualization, the FS 
within the upper contours are not plotted. 

Fig. 2 shows the Fermi surface used in our calcula- 
tion, and contours of = 23 meV. At this high en- 
ergy (about half of the maximum gap at antinode (it, 0), 
A^ ax = 44 meV), the 2-fold symmetry of the gap causes 
the contour pair approaching ends of the bananas to be 
more parallel near x-axis than the pairs near y-axis. The 
more parallel pair has better nesting effect locally (intra- 
or inter-contour) [13], which then further enhances the 
joint DOS . In view of Fig. 1(e) or 1(f), conspicuous in- 
terference peaks associated with wavevectors qi and q' x 
are clearly identified (see also Fig. 3). While, for exam- 
ple, no similar interference peak associated with q 2 (see 
Fig. 2) is seen. 

To illustrate how the qi and q' x peaks evolve as the 
change of bias energy, we plot D(q, u>) in Fig. 3, taking q 
along qa/n = (0,l)-(0,0)-(l,0) at bias energy from -5 to 
-23 meV (left: from bottom to top), and from +5 to +23 
meV (right: from bottom to top). Both frames show clear 
signatures of qi and q^ peaks which disperse as the bias 
energy changes. Their intensity increases as bias energy 



increases, to which qi disperses convergently to a fixed 
value near ±0.25(27r/a) y, while q' x disperses divergently. 
No peak develops at q 2 except some plateaus at interme- 
diate energies. The extensive, centered at zone center 
(0,0), quasi-lD like strong interference structure (paral- 
lel to k x axis, see Fig. 3 as well as Fig. 1(e)), comprised 
of q' x peaks discussed above, is superposed by the zero 
and small-q scattering connecting two k points within 
the same contour. It is the q^ peaks, which is enhanced 
by local nesting effect, gives rise to the ID interference 
structure parallel to k x axis. 
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FIG. 3. Scan of £>(q, ui) (arbitrary units) along qa/it — 
(0,l)-(0,0)-(l,0) at bias energy from -5 to -23 meV (left: from 
bottom to top), and from +5 to +23 meV (right: from bottom 
to top) at equal energy intervals. 

One particular interest on D(q, u) is the two qi peaks 
which have wavevectors near ±0.25(27r/a) y. These two 
peaks, if by themselves alone, would imply a ID parallel 
to x-axis stripe pattern in the real-space DOS (with a 
period of roughly 4a) — a pattern similar to the phe- 
nomenon predicted by "stripe" model [8]. 

To see how these stripes develop in real-space, in 
Fig. 4 we present the real-space LDOS modulation, tak- 
ing Fourier transform of D(q, u>). We consider only the 
uo = —23 meV and s — 0.2 case [Fig. 1(e)] and include 
the s = case for comparison. For s = 0.2 case, quasi- 
lD chains of beads of roughly 4a wide oriented approxi- 
mately along x axis are clearly seen at places away from 
the impurity site. The occurrence of these Friedel-typc 
oscillating stripes arises because of the strong, centered 
at (0, 0), ID interference structures. The latter also leads 
to the variation of intensities along y axis. In contrast, for 
the four-fold symmetry s — case, similar ID stripes also 
develop at places away from the impurity site. Clearly 
there exists no cross stripe pattern (checkerboard) any- 
where in the s = case, in contradiction to most theoret- 
ical predictions [5,6]. Recently Zhang, Hu, and Yu [14] 
considered the Zn or Ni impurity induced modulations of 
real-space LDOS for a d a .2_j / 2-wave superconductor which 
has a simple continuum band of circular FS. They [14] re- 



ported similar circularly one-dimensional (not cross) pat- 
terns. 

It should be emphasized that all the experimental q- 
space peaks are obtained upon Fourier transforming a 
single chosen piece of real space, which shows no sign of 
localized impurities. While in our approach, it is clear 
that the different q peaks are originated from different 
patches of space surrounding the impurity (see Fig. 4). 
To relate our calculation (s = case) to the experiments 
(checkerboard pattern), we need to have dilute impuri- 
ties, such that the decaying Friedel oscillations emanat- 
ing from individual impurities overlap crossly at some 
open spaces, and therefore Fourier transforming this sin- 
gle region gives all the q peaks. 




FIG. 4. Modulations of real-space LDOS for bias frequency 
u — — 23 meV. The left (right) frame is for a d a .2_ !/ 2+s-wave 
superconductor with s = (0.2) and both are presented for 
x,y — —40 to 40a. Single impurity is located at the center. 

A direct consequence of this picture is that when the 
impurity concentration is too high, the LDOS oscillations 
will be smoothed out, as the phases of the oscillations are 
pinned by individual impurities and therefore incoherent. 
This will be in contrast to the genuine stripe scenario [8] 
where the oscillation will be expected to be more promi- 
nent since the stripes, fluctuating in smooth real spaces, 
will now get pinned by the impurities. 

The case of multiple impurities for BSCCO is re- 
cently considered by Morr et al. [15] and Zhu et al. [16]. 
Very recently Anderson and Hedegard [17] also study the 
quantum interference between multiple impurities for a 
d x 2_ y 2-wave superconductor, whose pairing state could 
become d x 2_ y 2 +id xy or d x 2_ y 2 + is- wave at some doping 
level. Our dilute multi-impurity regime discussed above 
is fundamentally different from that any one else [15-17] 
has studied. Let the length scale between the impurities 
be do, and the distance away from the impurities, where 
we are interested, be d. Our interesting regime is then 
do ~ d > 20a (see Fig. 4), while the case other people 
have studied was do ~ d ~ a few a. 

Finally one does not expect any qualitatively new 
physics at the use of an unitary impurity, since it is known 
to affect only the LDOS locally close to the impurities, 
whereas we are now most interested in the regions away 
from the impurities. The essential point is that impurity 



scatters the quasiparticles, no matter in Born or unitary 
limit. 

In summary, motivated by recent ARPES and INS re- 
sults which are consistent with the picture that YBCO 
has superconducting gap of d x 2_ y 2+s symmetry, we pre- 
dict that quasi-one-dimensional Friedel-type stripes can 
be induced by defect in superconducting YBCO, observ- 
able by STM measurement. This phenomenon is most 
conspicuous at higher biases and associated with the lo- 
cal nesting effect for a Fermi liquid. Hopefully the pre- 
dicted modulations in either real or Fourier space can be 
tested by STM measurement soon. 

This work is supported by National Science Council of 
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